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Abstract 

A  new  semi-empirical  model  is  established  to  describe  the  cell  voltage  of  a  direct  methanol  fuel  cell  (DMFC)  as  a  function  of  current  density. 
The  model  equation  is  validated  experimental  data  over  a  wide  range  of  a  methanol  concentration  and  temperatures.  A  number  of  existing 
models  are  semi-empirical.  They,  however,  have  a  serious  mathematical  defect.  When  the  current  density,/,  becomes  zero,  the  equation  should 
reduce  to  the  open  circuit  voltage,  E0.  These  models,  however,  do  not  meet  the  mathematical  boundary  condition.  The  proposed  model  focuses 
on  very  unfavorable  conditions  for  the  cell  operation,  i.e.  low  methanol  solution  concentrations  and  relatively  low  cell  temperatures.  A  newly 
developed  semi-empirical  equation  with  reasonable  boundary  conditions  includes  the  methanol  crossover  effect  that  plays  a  major  role  in 
determining  the  cell  voltage  of  DMFC.  Also,  it  contains  methanol  activity  based  on  thermodynamic  functions  to  represent  methanol  crossover 
effect. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Characterization  of  the  fuel  cell  (FC)  frequently  uses 
large  and  complex  computer  models,  based  on  minute  de¬ 
tails  of  cell  component  design  (physical  dimensions,  mate¬ 
rials,  etc.)  along  with  chemical  and  physical  considerations 
(transport  phenomena,  electro-chemical  kinetics,  electrode 
kinetics,  etc.)  The  codes,  often  proprietary,  needed  in  the 
design  and  development  of  fuel  cells  are  cumbersome  and 
time-consuming  for  use  in  system  analysis  models.  Simpler 
approaches  are  normally  used  for  system  studies.  Another 
approach,  which  is  not  time  and  cost  efficient  would  be  to 
conduct  appropriate  tests  at  every  condition  expected  to  be 
analyzed  in  the  system.  Alternatively,  it  is  prudent  to  develop 
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correlations  based  on  the  thermodynamic  model,  which  de¬ 
scribe  cell  performance  according  to  operating  conditions 
such  as  temperature  and  pressure  [1]. 

In  the  development  of  model  equations  to  describe  the  per¬ 
formance  of  polymer  electrolyte  fuel  cells,  a  number  of  ap¬ 
proaches  using  empirical  models  have  been  attempted  [2-7]. 
In  many  cases  a  fairly  good  agreement  between  the  model 
and  experimental  data  is  achieved  by  adjusting  appropriate 
coefficient/parameters  in  the  model  equations.  Srinivasan  et 
al.  [2]  showed  that  it  is  possible  to  use  a  simple  equation 
to  describe  the  cell  voltage  vs.  current  density  behavior  for 
PEMFCs.  This  earliest  approach  tried  to  elucidate  the  behav¬ 
ior  of  such  a  complex  system,  and  was  solely  based  on  the 
electro-chemical  considerations,  which  formulated  an  empir¬ 
ical  equation  to  describe  the  polarization  curve.  This  empir¬ 
ical  equation  can  replicate  the  polarization  curve  reasonably 
well  at  low  and  intermediate  current  density,  but  failed  at 
high  current  density.  Kim  et  al.  [3]  improved  the  empirical 
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Nomenclature 

b  Tafel  slope 

E  cell  voltage 

Er  reversible  cell  potential 

Eq  open  circuit  voltage 

E*}  appropriate  open  circuit  voltage  (OCV) 

F  Faraday’s  constant 

J  current  density  (;) 

Jo  exchange  current  density  do) 

n  number  of  electrons  involved  in  a  reaction 

Re  ohmic  resistance  of  the  cell 

Greek  letters 

a  transfer  coefficient 

<5  methanol  crossover  effect  coefficient 

Ad  effective  coefficient  in  diffusion  overpotential 

for  mass  transport  limitation 


equation  for  better  replication  of  the  cell  performance  at  high 
current  density.  Squadrito  et  al.  [7]  reformed  the  Kim’s  equa¬ 
tion  with  addition  of  two  extra  terms  to  improve  the  prediction 
of  the  mass  transfer  related  resistance.  Argyropoulos  et  al.  [8] 
corrected  the  equation  that  coefficients  follow  specific  trends 
with  fuel  cell  operating  variables  and  allow  any  physically 
real  interpretation  of  the  model. 

Methanol  is  an  attractive  fuel  because  its  energy  density 
is  much  higher  than  that  of  hydrogen,  and  it  is  an  inexpen¬ 
sive  liquid  and  easy  to  handle,  store  and  transport.  How¬ 
ever,  in  practice,  DMFC  has  a  much  lower  open  circuit  volt¬ 
age  (OCV).  One  of  the  major  reasons  is  that  methanol  can 
cross  through  the  proton  exchange  membrane  (PEM),  such  as 
Nafion®,  to  reach  the  cathode  side  via  physical  diffusion  and 
electro-osmotic  drag  (by  protons).  Such  crossover  not  only 
results  in  a  waste  of  fuel,  but  also  lowers  the  cell  performance. 

The  effect  of  methanol  crossover  in  the  DMFC  has  attracted 
attention  worldwide  and  its  impact  on  cathode  operation  and 
system  efficiency. 

In  this  study,  we  develop  a  new  semi-empirical  model  for 
the  cell  voltage  as  a  function  of  current  density  response  of  the  > 

DMFC,  which  is  based  on  methanol  activity  with  thermody-  ^ 
namic  function,  diffusion  overpotential  to  represent  methanol 
crossover  contribution  and  mass  transport  limitation  which 
affects  significantly  the  performance  of  DMFC. 


2.  Model  development 

Srinivasan  et  al.  [2]  showed  that  it  is  possible  to  use  a 
simple  model  equation  to  describe  the  cell  voltage  ( E )  versus 
current  density  (j)  behavior  for  PEMFCs  in  the  activation  and 
ohmic  controlled  current  density  region: 

(1) 


Fig.  1 .  Vapor/liquid  equilibrium  experimental  data  for  water.  The  solid  lines 
are  calculated  from  this  work.  The  dotted  lines  are  calculated  from  the  PHSC 
model.  Open  squares  are  experimental  data  [11]. 


Current  Density(A  cm'2) 


Fig.  2.  Comparison  between  experimental  data  [8]  and  empirical  equation 
based  prediction  for  a  cell  operated  with  0.125  M  methanol  solution.  For  a 
cell  operated  with  each  0.125,  0.25  0.5  M  methanol  solution  supplied  at  a 
rate  of  1.12  cm3  min-1  with  air  fed  cathodes  pressurized  at  2  bar.  Catalysts 
consisted  of  35  wt.%  Pt-15  wt.%  Ru  (2  mg  cm-2  metal  loading)  and  mem¬ 
brane  type  is  Nafion  1 17  membrane  of  0.3  mm  thickness  (cell  temperatures: 
(♦)  343.15  K;  (A)  348.15  K;  (Q)  353.15  K;  (A)  358.15  K;  (0)  363.15  K). 


E=  E0-b  log  j  -  Rej 
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Current  density(A  cm'2) 

Fig.  3.  Comparison  between  experimental  data  [8]  and  empirical  equation 
based  prediction  for  a  cell  operated  with  0.25  M  methanol  solution.  For  a  cell 
operated  with  each  0. 125, 0.25, 0.5  M  methanol  solution  supplied  at  a  rate  of 
1.12  cm3  min- 1  with  air  fed  cathodes  pressurized  at  2  bar.  Catalysts  consisted 
of  35  wt.%  Pt-15  wt.%  Ru  (2  mg  cm-2  metal  loading)  and  membrane  type  is 
Nation  117  membrane  of  0.3  mm  thickness  (cell  temperatures:  (♦)  343. 15  K; 
(A)  348.15  K;  (O)  353.15  K;  (A)  358.15  K;  (0)  363.15  K). 

with 


E0=ET  +  b\ogj0  (2) 

where  Er  is  the  reversible  cell  potential,  b  the  Tafel  slope  for 
oxygen  reduction  and  R  the  ohmic  resistance  of  the  cell. 

Using  Eq.  (1),  with  the  appropriate  coefficients,  it  was 
shown  that  as  current  density  increased  the  predicted  cell 
potential  decreased  much  less  rapidly  than  observed  [2].  To 
increase  the  reliability  of  the  aforementioned  equation,  Kim 


Current  density(A  cm'2) 

Fig.  4.  Comparison  between  experimental  data  [8]  and  empirical  equation 
based  prediction  for  a  cell  operated  with  0.5  M  methanol  solution.  For  a  cell 
operated  with  each  0. 1 25, 0.25, 0.5  M  methanol  solution  supplied  at  a  rate  of 
1.12cm3  min~ 1  with  air  fed  cathodes  pressurized  at  2  bar.  Catalysts  consisted 
of  35  wt.%  Pt-15  wt.%  Ru  (2  mgcm~3  metal  loading)  and  membrane  type  is 
Nation  117  membrane  of  0.3  mm  thickness  (cell  temperatures:  (♦)  343.15  K; 
(A)  348.15  K;  (Q)  353.15  K;  (A)  358.15  K;  (0)  363.15  K). 


et  al.  [3]  suggested 

E  —  Eq  —  b  log  j  —  Rej  —  m  e”7  (3) 

where  m  and  n  are  the  parameters  that  account  for  the  “mass 
transport  overpotential”  as  a  function  of  current  density. 

Squadrito  et  al.  [7]  used  Eq.  (3)  as  a  starting  point  to  ana¬ 
lyze  the  different  contributions  to  the  mass  transport  limita¬ 
tion  and  produced  an  equation  in  the  form: 

E  =  E0  -  b  log  j  -  Rej  +  ajk  ln(  1  -  fij)  (4) 

where  a,  k  and  /3  are  the  parameters. 

The  term  Inf  1  —  fij)  introduces  a  limit  to  the  available  cur¬ 
rent  density.  For  k=  1,  a  has  the  same  dimension  as  Re  and 
can  be  interpreted  as  an  additional  resistance  term  due  to  the 
overall  mass  transport  limitation. 

Argyropoulos  et  al.  [8]  showed  the  applicability  of  Kim’s 
and  Squadrito’s  equations  for  predicting  voltage  response  of 
the  DMFC: 

Eceii  =  Eq  -  b  log  j  -  Rej  +  Cl  ln(  1  -  C2j)  (5) 


i.e.  in  the  Squadrito  equation  k=  0. 

Number  of  models  introduced  here  are  semi-empirical. 
All  models  are  based  on  Srinivasan  et  al.’s  model,  Eq.  (1). 
It,  however,  has  a  serious  mathematical  defect.  When  the 
current  density,  j,  becomes  zero,  the  equation  should  reduce 
to  the  voltage,  Eq.  These  models,  however,  do  not  meet  the 
mathematical  boundary  condition. 

We  propose  a  new  semi-empirical  cell  voltage  model: 

Eceii  =  E*  -  KRj  +  Xd  In  f  - — - - l\  (6) 

V  ,/luu  ) 

with 

Eq  =  Eq  =  <5<7MeOH  O) 


The  derivations  are  shown  in  Appendix  A.  We  employed 
methanol  activity  term  (<5«mcOh)  to  take  into  account  the 
methanol  crossover  effect  account  [9]  and  S  is  an  effec¬ 
tive  coefficient.  To  describe  methanol  activity  of  DMFC,  we 
employed  an  equation  of  state  (EOS)  based  on  the  mod¬ 
ified  perturbed  hard-sphere-chain  (PHSC)  model.  It  con¬ 
siders  a  new  perturbation  equation  that  is  obtained  from 
the  generalized  Lennard-Jones  potential  function  based  on 
statistical-mechanical  relationship  [10]: 


In  a  = 
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Table  1 

Calculated  values  for  the  this  semi-empirical  equation  coefficients  for  DMFCs  operated  with  methanol  solution 


Cell  temperature  (K) 

MeOH  solution  concentration  (M) 

Calculated  values 

Eq  (V) 

MV) 

arPCT1) 

M(V) 

Limiting  current 

363.15 

0.125 

1.3175 

2.6557 

6.8525 

0.07097 

0.0421 

358.15 

0.125 

1.2876 

2.6328 

7.5806 

0.06451 

0.0337 

353.15 

0.125 

1.2708 

2.6103 

6.4519 

0.07327 

0.0282 

348.15 

0.125 

1.2579 

2.5882 

5.5692 

0.0921 1 

0.0230 

343.15 

0.125 

1.2460 

2.5665 

7.6569 

0.07663 

0.0191 

363.15 

0.25 

1.2238 

2.3657 

1.9491 

0.06022 

0.0939 

358.15 

0.25 

1.2136 

2.3244 

2.9761 

0.02846 

0.0870 

353.15 

0.25 

1.2021 

2.2832 

3.6132 

0.04775 

0.0731 

348.15 

0.25 

1.1934 

2.2412 

3.1202 

0.06307 

0.0591 

343.15 

0.25 

1.1538 

2.2008 

1.8737 

0.06749 

0.0482 

363.15 

0.5 

1.1886 

2.1041 

1.8067 

0.02090 

0.2001 

358.15 

0.5 

1.1662 

2.0832 

1.7101 

0.02849 

0.1702 

353.15 

0.5 

1.1650 

2.0598 

1.7872 

0.03399 

0.1581 

348.15 

0.5 

1.1423 

2.0357 

1.8624 

0.03591 

0.1278 

343.15 

0.5 

1.1278 

2.0125 

2.0314 

0.03451 

0.1152 

The  first  term  of  Eq.  (6),  Eq  is  the  voltage  without  methanol 
crossover.  In  Eq.  (7),  Eq  is  defined  as  an  appropriate  open 
circuit  voltage  (OCV).  In  this  study,  when  the  current  density 
becomes  zero,  all  terms  should  be  zero  except  £q  term.  In 
Eqs.  (1),  (3)-(5),  the  term,  log j,  however,  goes  to  infinity. 
It  is  not  possible  to  have  infinity  at  j  =  0.  The  third  term,  atr/' 
presents  the  ohmic  resistance  and  electrode-electrolyte  over¬ 
potential  of  the  cell.  To  take  into  account  rapid  voltage  drop 
at  high  current  density,  we  employ  diffusion  overpotential  in 
last  term. 

The  above  equation  describes  the  methanol  crossover  con¬ 
tribution  and  diffusion  contribution  for  mass  transfer  limita¬ 
tion  region.  In  general,  the  methanol  crossover  effect  plays  a 
major  role  in  determining  the  open  circuit  voltage  of  DMFC. 
S  and  aj  represent  the  methanol  crossover  effect  coefficient 
and  effective  coefficient  in  diffusion  overpotential  for  mass 
transport  limitation,  respectively.  Also  value  of  jyim  is  value 
of  j  when  i]  — »  —  oo. 

Confirmed  in  Eq.  (6),  as  the  value  of  the  current  den¬ 
sity  approaches  zero,  the  cell  voltage  becomes  the  appro¬ 
priate  open  circuit  voltage  (OCV),  Eq  containing  methanol 
crossover. 


3.  Result  and  discussion 

We  propose  a  new  semi-empirical  model  to  predict  the 
cell  voltage  as  a  function  of  a  current  density  for  a  liquid 
feed  direct  methanol  fuel  cell  (DMFC).  The  performance  of 
the  DMFC,  at  higher  current  densities  has  been  shown  to  be 
limited  by  a  mass  transport  process,  which  mainly  occurs,  by 
diffusion. 

In  Fig.  1,  the  solid  and  dotted  lines  are  calculated  from 
the  proposed  model  and  MPHSC  equation,  respectively  [11]. 
In  the  entire  density  region,  MPHSC  EOS  shows  slight  devi¬ 
ations  from  the  experimental  data  [11],  however  our  model 
shows  good  agreement  with  the  data.  The  obtained  charac- 
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Fig.  5.  Experimental  overvoltage  values  of  only  methanol  crossover  [9] 
(temperature:  348  K,  pressure:  2.0  bar  air,  methanol  concentration:  0.125  M). 


Fig.  6.  Experimental  overvoltage  values  of  only  methanol  crossover  [9] 
(temperature:  348  K,  pressure:  2.0  bar  air,  methanol  concentration:  0.5  M). 
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teristic  parameters  for  water  are  rMeOH  =  2.9,  crMeOH  =  2.95, 
£MeOH/^=  550.2. 

In  Figs.  2-4,  calculated  cell  voltages  are  presented  for 
three  different  aqueous  methanol  solution  concentrations 
(0.125,  0.25  and  0.5  M)  at  a  various  range  of  cell  operat¬ 
ing  temperatures.  Model  parameters  are  listed  in  Table  1. 
The  model  equation  proposed  in  this  study  is  able  to  ex¬ 
press  the  cell  voltage  in  the  entire  current  density  region.  The 
existing  models,  however,  do  not  satisfy  the  mathematical 
boundary  condition,  that  is,  when  the  current  density  goes  to 
zero,  the  open  circuit  voltage  (OCV)  diverges  to  the  infin¬ 
ity.  However  our  calculated  Eq  values  are  reasonable  with 
temperature  as  shown  in  Table  1.  In  practical,  the  value  of 
OCV  decreases  with  increasing  methanol  concentration  [9]. 
As  shown  in  these  figures,  the  cell  voltage  decreases  with  in¬ 
creasing  current  density.  Especially,  at  the  high  current  den¬ 
sity  region,  the  cell  voltage  decreases  radically.  It  is  because 
the  diffusion  affects  the  performance  of  total  cell  voltage. 
Thus,  we  employ  these  two  contributions  in  the  proposed 
model. 

Undoubtedly,  the  methanol  crossover  is  not  the  only  is¬ 
sue  in  DMFC.  One  must  also  consider  aspects  such  as  the 
pressure  effect  and  fuel  flow  rate.  In  fuel  cells,  membranes 
operate  under  the  significant  mechanical  pressure  that  acts 
in  opposition  to  the  osmotic  pressure  and  fuel  flow  rate  is 
involved  in  a  methanol  crossover. 


4.  Conclusion 

We  have  established  a  new  semi-empirical  model  to  de¬ 
scribe  the  cell  voltage  of  a  direct  methanol  fuel  cell  (DMFC) 
as  a  function  of  current  density.  The  model  presented  here 
includes  an  accurate  open  circuit  voltage  (OCV),  methanol 
crossover  from  methanol  activity  and  diffusion  effect  that 
takes  into  account  rapid  voltage  drop  at  high  current  density. 
Its  advantage  follows  from  its  simplicity;  a  simple  algebraic 
form  with  a  few  adjustable  model  parameters  appears  to  be 
suitable  for  representing  the  appropriate  open  circuit  voltage 
for  the  given  systems. 


where  E0 cv  is  the  open  circuit  voltage,  FmethanoLcrossover 
the  methanol  crossover  effect,  i7eiectrode-eiectrolyte  the 
electrode-electrolyte  overpotential,  ??ohmic  .resistance  the  ohmic 
resistance,  ^diffusion  the  diffusion  overpotential. 

A.  1.  Methanol  crossover 


The  methanol  crossover  effect  is  expressed  by  methanol 
activity  equation. 

As  shown  in  Figs.  5  and  6,  the  methanol  crossover  can  be 
represented  by: 

methanol  crossover  —  <5«MeOH  (A2) 


A.2.  Electrode-electrolyte  overpotential 


Butler- Volmer  equation  is  as  follows: 


l  =  IQ 


exp 


n(\  —  a)Fi]\ 


RT 


'  —naFn\ 

expi^H 


(A3) 


where  i  is  the  current  density  (it  has  same  meaning  of  jj,  io 
the  exchange  current  density  (jo),  a  the  transfer  coefficient, 
n  the  number  of  electrons  involved  in  a  reaction  and  F  the 
Faraday’s  constant. 

Taylor  expansion  gives: 


J  —  Jo 


n(\-  a)F  ( n(  1  -  a)F\~  2 

1  +  — ™ — '7  +  — ^ —  r  ■ 


RT 


—naF 

—  <  1  H - r\ 

I  RT  1 


RT 

2 


naF\  9 

)  r 


RT 


=  Jo 


n(\  —  a)F  —naF  j 

RT~  )  11 


RT 


n(  1  —  a)F 
RT 


- naF 
RT 


(A4) 


We  take  first  order  approximation  due  to  77  <§;  1 : 
RT 

^electrode— electrolyte  —  77rJ 
nFjo 


(A5) 
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Appendix  A 

The  total  cell  voltage  is  given: 

Cecil  —  FqcV  ^methanol  crossover  ^electrode— electrolyte 

(Al) 


Ohmic  resistance  is: 

Pohmic  jesistance  —  R  c  j  (A6) 

where  Re  is  the  ohmic  resistance  of  the  solid  polymer  elec¬ 
trolyte. 

A.  4.  Diffusion  overpotential 

Nernst  diffusion  layer  equation  is: 
n  RT 

E  —  E°  -) - In  ao 

nF 


7?ohmic  .resistance  ^diffusion 


(A7) 
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Above  equation  can  represent  current  density: 


By  combining  Eqs.  (Al),  (A2),  (A5),  (A6)  and  (A13): 


n  RT 

Ei  =  E°  H - In  a 

nF 


(A8) 


Diffusion  overpotential  is  difference  of  two  overpotential: 


£cell  =  Eq~  KR  j  +  kd  In 


with 


(A14) 


n  =  Ej  —  E  =  — — 


RT  a 

E  =  - In— 

nF  a  o 


The  activity  can  be  represented  concentration  ratio: 
'  n  Fvj  \ 


C  r0 
—  =  —exp 
C0  r 


RT  ) 


(A9) 


(A10) 


Current  density  can  be  expressed  combining  Fick’s  law  about 
diffusion  and  Faraday’s  law: 

C0-C 


—  =  —D- 

nF 


l 


(All) 


Combining  (A10)  and  (Al  1),  current  density  is: 


J 


n  FDCq 

1 

n  FDCq 

7 


r  o 


nFrj 


1 - exp  - 

'  RT 


nFrj 

‘"“Pi  Hr 


=  jlim  1  -  exp 


nFr] 

RT 


(A12) 


where  as  jiim  is  limiting  current  density,  it  is  value  that 
1]  — >  — oo. 

Therefore,  diffusion  overpotential  can  be  derived  as  fol¬ 
lows: 

^diffusion  =  “In  (  7  '  )  =  koln  f  7  -  J  (A13) 

nF  \  J lim  J  \  Jlim  J 


Eq  =  Eo  —  <5flMeOH 


(A15) 


KR  j 


J 


(Alb) 


References 

[1]  K.  Kordesch,  G.  Simader,  Fuel  Cells  and  Their  Applications,  VCH 
Verlagesellschaft  mbH,  Weinheim,  1996. 

[2]  S.  Srinivasan,  E.A.  Ticianelli,  C.R.  Derouin,  A.  Redondo,  J.  Power 
Sources  22  (1988)  359. 

[3]  J.  Kim,  S.  Lee,  S.  Srinivasan,  C.E.  Chamberlin,  J.  Electrochem.  Soc. 
142  (1995)  2670. 

[4]  J.H.  Lee,  T.R.  Lalk,  J.  Power  Sources  73  (1988)  229. 

[5]  J.H.  Lee,  T.R.  Lalk,  A.J.  Appleby,  J.  Power  Sources  70  (1988) 
258. 

[6]  J.H.  Lee,  T.R.  Lalk,  A  Modeling  Technique  for  Fuel  Cell 
Stack  Systems,  ASME  Dynamics  Systems  and  Control  Division, 
1996. 

[7]  G.  Squadrito,  G.  Maggio,  E.  Passalacqua,  F.  Lufrano,  A.  Patti,  J. 
Appl.  Electrochem.  29  (1999)  1449. 

[8]  P.  Argyropoulos,  K.  Scott,  A.K.  Shukla,  C.  Jackson,  Fuel  Cells  2 
(2002)  78. 

[9]  K.  Scott,  W.M.  Taama,  P.  Argyropoulos,  K.  Sundmacher,  J.  Power 
Sources  83  (1999)  204. 

[10]  J.Y.  Seong,  Y.C.  Bae,  submitted  for  publication. 

[11]  J.M.  Smith,  H.C.  Van  Ness,  M.M.  Abbott,  Chemical  Engineering 
Thermodynamics,  McGraw-Hill,  2002. 


